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Abstract. Several observational analyses suggest that 1. Introduction 



matter is spatially structured at a scale of ^lSS ~ 
130 h^^ Mpc at low redshifts. This peak in the power spec- 



trum provides a standard ruler in comoving space which parent et al. 1986; Geller & Huchra 



It has been known for more than a decade (e.g. ie Lap- 
that the spa- 



can be used to compare the local geometry at high and low 
redshifts, thereby constraining the curvature parameters. 

It is shown here that this power spectrum peak is 
present in the observed quasar distribution at z '--^ 2: quali- 
tatively, via wedge diagrams which clearly show a void-like 
structure, and quantitatively, via one-dimensional Fourier 
analysis of the quasars' tangential distribution. The sam- 
ple studied here contains 812 quasars. 

Th e method produces strong constraints (68% confi - 



tial distribution of extragalactic objects is structured at 
length scales about an order of magnitude greater than 
the ro « 5 Mpc scale which characterises galaxy clus- 
tering via the two-point auto-correlation function. Obser- 
vational analyses from several different data sets suggest 
that this is due to a characteristic length scale at ^Lgg ~ 
130 h^^ Mpc, or in other words, that there is a maximum 
in t he power spectrum at fc = 27r/£T-c;s ^ . 05/tMpc~ "^ 



(e.g. 



Broadhurst et al. 1990; Broadhurst 1999; Broadhurst 



dence l imits) on the dcnsitjr parameter Qq and weaker con - 



Jaffc 1999; ia Costa 1992 



straintt on the cosmological constant Aq^ which can b 



Efstathiou 1993, 1994 



pressed 



Indepe ndently of Aq (in the range Aq £ [O;!]); the con - 



by the relation np- (o72d±OT5)-H(o7lO±0^^ ^\ llggzii 



ia Costa et al. 1993 ; Baugh fc 



Gaztanaga fc Baugh 1998; Einastc 



Deng, Deng fc Xia 1996| or Guzzo 



1999 for a recent review) 



straint is 0.1 < fto < 0.45. Constraints if the cosmological 
constant is zero or if Aq = 1 — ilo arc Hq = 0.24lQ "g and 
flo = 0.30 ± 0.15 respectively. 

The power spectrum peak method is independent from 
the supernovae Type la method by choice of astrophysical 
object, by redshift range, and by use of a standard ruler 
instead of a standard candle. Combination of the two re- 
sults yields f^o = (0.30 ± 0.11) + (0.57 ± 0.11)(Ao - 0.7), 
0.55 < Aq < 0.95, (68% confidence limits) without assum- 
ing that Xq = 1 — Qq. This strongly supports the possi- 
bility that the observable universe satisfies a nearly flat, 
perturbed Friedmann-Lemaitre- Robertson- Walker model, 
independently of any cosmic microwave background obser- 
vations. 



Key words: cosmology: observations — cosmology: the- 
ory — distance scale — quasars: general — large-scale 
structure of Universe — reference systems 



It has already been suggested that this scale could be 
used as a standard ruler which could be compared at low 
and high redshifts in order to constrain the curvature pa- 
rameters, rig (the density parameter) and Aq (the dimen- 
sionless cosmological constant). In the redshift direction, 
at least two analyses have been carried out based on this 
idea: one analysis of a quasar catalogue and one analysis 
of Lyman break galaxies. 

Deng, Xiaoyang & Fang (1994) implicitly used the 
large scale structure scale as a curvature constraint in 
quasar data, under the assumption that Aq = 0, and found 
that fio ~ 0.4. Broadhurst & Jaffe used the radial (red- 
shift) distribution of Lyma n brea k galaxies at 2; ~ 3 (ta- 
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ble 1, iGiavahsco et al. 1998] ; fig. 2, [Adelberger et al. 1998|) . 
They found a correlation scale of Az « 0.22 ± 0.02, and 
inferred a relation 3.2 51o ~ ^0 ~ 0.7. 

The purpose of this paper is to (i) emphasise that the 
principle can be applied to a class of bright objects eas- 
ily found at super-unity redshifts: quasars; (ii) show that 
redshift selection effects can be minimised by using the 
tangential distribution instead of the radial distribution; 
and (iii) show pictorially (i.e. qualitatively) that quasars 
do indeed trace large scale structure at z ^ 2. 



2 



B. F. Roukema & G. A. Mamon: Curvature Parameters from QSO Tangential LSS 



(i) Because quasars are much brighter than Lyman 
break galaxies, they offer a potentially much more rapid 
method of obtaining high precision estimates of the cur- 
vature parameters than the latter. Both classes of objects 



have the advantage relative to supernovae type la (Perl- 
mutter ^F^1~T99^; Riess et al. 1998) of being at super- 
unity redshifts, so that the dependence on the curvature 
parameters is strong. 

(ii) In order to avoid the well-known selection effects 
in the redshift distribution of quasars (which could also in 
principle affect the redshift distribution of Lyman break 
galaxies), the tangential distance distribution of quasars 
is investigated. 

(iii) However, two-dimensional wedge diagrams are 
also plotted to show clear qualitative evidence of the trac- 
ing of large scale structure by the quasar distribution. 

Details of the method and selection of homogeneous 
quasar samples are described in Sect. ||. Results are pre- 
sented in Sect. 3, and discussion and conclusions are pre- 
sented in Sect. 4. 

A perturbed Friedmann-Lemaitre-Robertson- Walker 
cosmological model is assumed here. The context in which 
quasars could reasonably be expected to form a tracer 
population of large scale structure may be any model in 
which quasars form in galaxy centres or in which galax- 
ies form around quasars. The expected short life times of 
quasars should not prevent them from forming a tracer 
population, though they might form a biased population 
which could either weaken or strengthen the amplitude of 
the signal. 

The Hubble constant is parametrised here as ft. = 
TJo/lOOkm s~^ Mpc~^. Comoving coordinates are used 
throughout [i.e. 'proper distances', eq. (14.2.21), Weinberg 
(1972|2 equivalent to 'conformal time' if c = 1]. Values of 
the density parameter, fioj and the dimensionless cosmo- 
logical constant, Aq, are indicated where used. 




Fi g. 1. Redshift distribution of quasars in the compilation 
of [Vcron-Cctty fc Vcron (1998| ). 
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2. Method 

2.1. Choice of catalogue and sky regions 

There are now around 10^ quasars which have publicly 



available redshifts and celestial positions (Veron-Cetty & 
Veron 1998[ ). What is the optimal way to search for the 
power spectrum peak among these data? 

Figs 1^ and |2| show the redshift and sky distributions 
of these quasars. Depending on the redshift range of any 
quasar sample, a few hundred Mpc typically corre- 
spond to Az ~ 10""'^ in redshift to within an order of 
magnitude, and at z ^ 2, where the peak in the redshift 
distribution lies, to A9 ^ 1°. 

There is clearly structure in the combined redshift his- 
togram for the full sky at around this scale. As discussed 



in detail by Scott (1991), the ratios in wavelength of im- 



portant emission lines which contribute to the chance of 
detecting a quasar (Lya, C IV, C III, Mg II) correspond 



Fig. 2. Sky distribution of quasars in the compilation of 
Veron-Cetty &: Veron (1998|) . 



to intervals AlnA w 0.2 and clearly contribute to the ob- 
vious peaks in the distribution. 

The alternative to searching for structure in the red- 
shift direction is to search in the tangential direction. Most 
of the deep surveys visible in Fig. || are based on photo- 
metric selection from photographic plates of size roughly 
6° X 6°, in particular from objective prism surveys. This 
is moderately larger than the scale of interest. 

In order to minimise possible systematic effects due 
to obscuration by dust and modification of the sky back- 
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Fig. 3. Redshift distri bution of quasars with redshifts 
from [ovino et al. (199^ ). 



Table 1. Definition of the two tangentially long, densely 
observed, homogeneous, subsamples of the lovino et al 



(1996}_ quasar survey listed in the Veron-Cetty fc Veron 
(1 998 catalogue, near the south galactic pole (SGP). The 
right ascension (first row) and declination (second row) 
subsamples are defined by J2000 limits {ai < a < a2, 
Si < S < 32)- For the purpose of Poisson simulations, 
these are subdivided in right ascension at {a', a") and in 
declination at {d',S"), in order to allow for the possibility 
of different magnitude zero points or different magnitude 
cutoffs in the different plates. The number of objects N in 
each subsamplc is indicated. The total number of physi- 
cally distinct quasars in the two subsamples is A^tot = 812. 



cei a.1 &\ &2 01 a" 


5' 


5" N 


'Right ascension (a) subsample' 

Qh42™ l'^59™ -42.0 -37.5 I'^OT™ l''33'" 

'Declination (5) subsample' 

Qh^^m ^^00™ -42.0 -28.0 


-37.5 


604 

-32.5 373 



ground by bright stars, surveys near the South or North 
galactic poles would be best. 

Near the South galactic pole (SGP), several regions 
have been observed contiguously, with (at least to the eye) 
a reasonable homogeneity across the different plates in the 
region. 

More o bjectively, the largest single h omogeneous sub- 
set of the Veron-Cetty fc Veron (1998) catalogue in the 
SGP region is that of lovino, Clowes fc Shaver (1996 ), who 




SO'" IS"" 
a (J2000) 



Fig. 4. Sky distri bution of quasars with redshifts from 
lovino et al. (199*^ ), which lie in the redshift range 1.8 < 
z < 2.4. Thick lines show the outlines of the right ascen- 
sion and declination samples chosen for one-dimensional 
Fourier analysis. Thin vertical and horizontal lines show 
divisions in right ascension and declination in the right 
ascension and declination subsamples respectively, which 
simulated distributions are Poisson distributed. See Ta- 
ble |l| for numerical values of these limits. 

used an 'automatic quasar detection' method, i.e. applied 
a computer algorithm to the digitised images of objective 
prism plates taken on the UK Schmidt Telescope at Siding 
Spring, Australia. 

This is the catalogue chosen for analysis here. The 
overall redshift distribution of this catalogue is shown in 
Fig. ^, and the sky distribution of the 1.8 < z < 2.4 com- 
ponent is shown in Fig. |^. 

Wedge diagrams of the catalogue, within the limits of 
Table |l| (see Sect. are shown in Fig. |. Note that 
the right ascension and declination subsamples are not 
entirely independent sets of quasars (see Fig. U), though 
since only one angular coordinate is used in each analy- 
sis, they are very close to being effectively independent. 
Although quantitative analysis in these planes would be 
difficult to carry out due to the obvious redshift selec- 
tion effects, it is qualitatively clear that a void-like struc- 
ture is present, at a scale near ^Lgg ~ 130 Mpc for 
(r^o — 0.3, Ao = 0.7 ). In order to avoid the redshift selec - 
tion effects (see also Hartwick fc Schade 199C| ; pcott 1991 ), 
the analysis here is restricted to Fourier space. 



2.2. One- dimensional Fourier Analysis 

The study of structure at ^lss is made by one- 
dimensional Fourier analysis of two subsamples of the 



4 



B. F. Roukema & G. A. Mamon: Curvature Parameters from QSO Tangential LSS 




X (h"'Mpc) 




X (h"'Mpc) 



Fig. 5. Wedge diagrams for the for the right ascen- 
sion (upper panel) and declination (lower panel) subsam- 
ples, for (r^o = 0.3, Ao = 0.7). Since curvature is zero, 
rectlinear coordinates are defined: x 



y ^ d 



pm 



sm f " - 



where 0* = 



^pm 



COS 



62) /2 + 7r/2 and 



9, 61,62 and dpm are defined in Sect. |2.2| . Voids consistent 
with ^LSS = 130 ± 10 Mpc and not due to selection 
effects are clearly visible. 



lovino et al. (1996) sample [i.e. of the subset of the Vcron- 



Cetty ( 'z Veron catalogue for which rcdshifts are obtained 
from lovino et al. (1996| )] which have maximum tangen- 
tial survey length. That is, one subsample combining three 
plates in the right ascension direction and one subsam- 
ple combining three plates in the declination direction are 
chosen. 

The right ascension and declination boundaries are 
chosen conservatively, i.e. to within at most 2.5 great cir- 
cle degrees from the plate centres, and are shifted even 
closer in where it looks like there may be incompleteness 
close to the boundary. Although what appears to be a lack 
of quasars near a few plate boundaries might in fact be 
due to real voids, it is preferable to risk losing some real 
signal rather than risk including some noise. 



The angular limits chosen are indicated in Table ^ 
The redshift range used iszi = 1.8 < z < Z2 — 2.4, 
which includes most of the catalogue, and is small enough 
to superimpose only a few 'units' of large scale structure, 
i.e. the signal should not be significantly reduced by the 
superimposition of structures which are out of phase. 

In each of the two subsamples, the angular positions 
in the long direction (right ascension and declination re- 
spectively) are converted to comoving tangential lengths 
dj_ by 



dj_{z, t 



61) dpm{z) 

Rc smh[d(z) / Rc], kq < 

7i) ■( d{z), Ko = 

Rcsm[d{z)/Rc], kq > 0. 



(1) 



where the angular position 6 ~ a cos d or 6 ^ S and the 
survey limit is 61 = ai cos 61 or 61 = Si respectively (in 
radians), dpm(-z) is the proper motion distance. 



diz) 



c 

i5o 



da 



/(i+z) ay/Qo/a - Ko + Xqc 



(2) 



is the proper distance [eq. (14.2.21), Weinberg (1972)] 



Kq = Hq + Xq - 1 



(3) 



500 R 



is the (dimensionless) curvature of the observational Uni 
verse and 

c 1 



(4) 



is its curvature radius. 

The maximum tangential length considered is 
d±{zi,d2), where zi = 1.8 is the low redshift limit of 
the sample. This is the minimum tangential length cor- 
responding to 62 — 61 over the range in z and the {flo, Ao) 
domain considered in this paper. For negative or zero cur- 
vature, d± is always an increasing function of z. For pos- 
itive curvature and flo < 1.1, Aq < 1, the d/Rc ~ tt/2 
point (halfway to the antipode) occurs at z >, 8, so the do- 
main of decreasing d± is not reached here. Thus, d± is a 
strictly increasing function of z in the domain of interest 
of this paper and d±{zi) provides the minimum tangential 
length. 

Note that the choice of this cutoff throws away a small 
amount of data (e.g. for which z ^ Z2,d ~ 62), but en- 
sures that the one-dimensional number density distribu- 
tion dN/dd± is a uniform projection of a (large) subset 
of the two-dimensional distribution d^iV/ddj_ddpm. In- 
clusion of the small amount of lost data would create a 
nonuniform number density projection and would intro- 
duce non-physical power to the Fourier transform, so is 
not attempted here. 

For values of (fioiAo) in the range 0.0 < ^0 < 
1.1,-0.1 < Ao < 1.0, the list of positions < dj_ < 
d±{zi, 62) is binned into 1024 bins and fast Fourier trans- 
formed to a function /(z^jilo, Aq). The results below are 
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found to change slightly but insignificantly if fewer bins, 
e.g. 128 bins, are used. The contours in the (flo, Aq) plane 
are more noisy with fewer bins. 

Two null hypotheses are considered here for the 
Fourier spectrum for each pair (fiojAo): the possibility 
that large scale structure is undetectable in the two quasar 
subsamples, and the possibility that the best estimate of 
the scale of large scale structure is ^lSS ^ ^-^LSS- 

2.3. Null Hypothesis Tii: No LSS peak is detectable 

Informally, Hi{flo, Xq) is the hypothesis that there is no 
peak in the tangential distributions of the quasars due to 
large scale structure. 

More precisely, Hi(rio,Ao) is the hypothesis that: 

(i) the pair (fioj^'^o) is correct; 

(ii) the value of the Fourier transform / at the large scale 
structure frequency vq = 1/LlsS' i-^- /(I'ol^^o, Aq), is 
not significantly higher than that expected from the 
distribution of the same statistic evaluated for Poisson 
distributions in 9, for fixed values of i^o, ^^Oi Aq, for both 
subsamples; and 

(iii) z^max, defined as the local maximum in f{i'\Qo,Xo) 
at the greatest value of satisfying v < 0.01 /iMpc~^, 
is not significantly higher than that expected from the 
distribution of the same statistic evaluated for Poisson 
distributions in 9, for fixed values of rio,Ao, for both 
subsamples. 

The Poisson distributions are pseudo-random sam- 
plings of uniform distributions in 9 within each of the 
three subdivisions of the sample as defined as in Table |l|. 
The number of simulations calculated is 30. A Gaussian 
smoothing of standard deviation two bin widths is applied 
to the Fourier transform /{vIHq, Aq) before searching for 
the local maximum. 

The purpose of criterion (ii) is that if a peak is present 
at the scale expected, then this may contribute to rejecting 
Hi by the presence of a strong peak. However, because 
the overdensity may not be very high, this may not be 
sufficient in itself to reject Hi. 

An independent and possibly more sensitive test is 
(iii): is the best estimate of the frequency of a peak signif- 
icantly different from that for Poisson distributions, in- 
dependently of any criterion on the absolute height of 
the peak? Likely values of (fmax) for Poisson simulations 
are around 0.005 /i Mpc~^, though this depends on the 
smoothness or roughness of / for the Poisson simulations. 
Since l/iLgg w 0.0077 /i Mpc~^, then, as long as the scat- 
ter in Vmax for the simulations is small enough, criterion 
(iii) may enable rejection of Hi for pairs (r2o,Ao) which 
correctly describe the observational Universe, because in 
that case (under the principle assumed for this paper) , the 
peak is expected to occur at a special frequency rather 
than at an arbitrary frequency. 



2.4- Null Hypothesis 7^2 • The best estimate of the 
frequency of an LSS peak is at 1/LLgg 

The more interesting hypothesis is H2i^o, Xq), the hy- 
pothesis that the best estimate of a peak in the tangential 
distributions of the quasars due to large scale structure is 
at i^Lss ^ ^^LSS' independently of whether the peak is 
significant or not. 

7^2 (f^O: Ao) is quantified as follows: 

(i) the pair (r2o,Ao) is correct; 

(ii) t'max (defined as for Hi) is consistent with fmax = 
1 / ^Lss where the 1 a uncertainty in the external es- 
timate of .Z^Lgg is Ai^Lgg, and Ai^rnax, the la un- 
certainty in estimating z^max from the present data is 
obtained robustly by bootstraps (e.g. Barrow, Sonoda 



fc Bhavsar 1984) 



The estimate AL^gg = 10 h^^ Mpc is adopted here. 

The bootstrap method (e.g. Barrow et al. 1984) for a 
catalogue of N objects, is for N objects to be randomly 
drawn from the same sample, allowing multiple sampling 
of single objects. The statistical uncertainties in the prop- 
erties of interest are then estimated by running several 
such bootstrap simulations, which are considered as inde- 
pendent experiments. This provides an upper estimate to 
the uncertainty. The number of bootstraps used here is 
30. 

The bootstrap 1 a uncertainty and AL^gg are as- 
sumed to be independent and to arise from Gaussian dis- 
tributions, so are combined in quadrature. 

Although H2 could in principle be consistent with the 
data for (f2o: Ao) pairs which are also consistent with Hi, 
the (f2o,Ao) pairs for which Hi is rejected and H2 is not 
rejected are obviously of most interest. 

3. Results 

3.1. Hi: Can the absence of a peak be rejected? 

Confidence levels 1- Pa[f{i^o)], l--Pa(t'max), l-Psifit^o]] 
and 1 — Ps (i^max), for rejecting Hi are shown in Figs g 
and for the right ascension and declination subsamples 
respectively. 

These are defined by the probability of the observa- 
tional results given the hypothesis: 



P{t)=Pganssit,t,At) = J 



t-t| /At 



/27r 



du, 



ierfcf^^ 



(5) 



where t is the parameter studied (either f{i'a) or i^max), 
and t and Ai are the mean value of t and the standard 
deviation of t obtained from the Poisson simulations. 

Since the question of interest is to find (f2o, Aq) pairs 
for which there is, ideally, an excess of power at ^Lgg 
and a frequency higher than that for Poisson distributions 
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Fig. 6. Confidence levels for the different criteria for try- 
ing to reject TLi, for the right ascension subsample. The 
upper panel is for fii^o), the lower panel for i^max- Rejec- 
tion of the null hypothesis Hi at > cr, > 1 cr, > 2 cr and 
> 3(7 confidence levels, i.e. 1 - P > 50%, 1 - P > 84%, 
1 - P > 98% and 1 - P > 99.9%, is shown by contoured 
regions with light shading, light cross-hatched, medium 
cross-hatched and heavy cross-hatched shading, respec- 
tively. (In a few of this set of figures, a few contours at 
< Ofj and < —1(7 are also shown. These are not useful 
for null hypothesis rejection.) Lines indicating Ao = and 
Ao = 1 — r^o are shown for f2o < 1 as a guide to the eye. 



(given t he li mits for searching for a local maximum defined 
in Sect. 



2.3), the probability above is defined to be one- 



sided. 

The upper panel of Fig ^ shows that if there is a peak 
at ^LSS right ascension subsample, then it is not 

strong enough to significantly reject the null hypothesis 
Tiiiflo, Xq) of the non-existence of a peak for any pair 



Fig. 7. Confidence levels for rejecting TCi, for the decli- 
nation subsample. The upper panel is for /(I'o), the lower 
panel for t'max- Shading is as for Fig. |. 



{i^Q,Xo). Note, of course, that non-rejection of Hi does 
not imply that Hi is correct. It just implies (states) that 
Hi is not rejected. 

However, the lower panel of Fig ^ shows that the best 
estimate of the frequency of a peak, independently of its 
significance, is rejected at the 1 — P > 84% level for a 
large (though noisy) band in the {flo, Aq) plane, for flo ~ 
0.4 ±0.1. 

Both the /{vq) amplitude test and the z/max test for 
the declination subsample (Fig |^) independently confirm, 
to ~ ±0.2 precision in and Ao, the region of the (^o, Ao) 
plane which enables la rejection oi Hi{^lo, Xq). A small 
band within the latter, for the Vmax test, enables rejection 
to 2 a, i.e. 1 - P > 98%. 

As noted above, the fact that the /(z^o) test for the 
right ascension sample does not reject Hi is not a problem 
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Fig. 8. Combined confidence levels for rejecting 7i i , using 
Eq. to combine the results shown in Fig. || and Fig. ^. 
Shading is as for Fig. |. 

for alternative hypotheses to Hi (in particular for 0.2)- 
The failure of one test to reject a model does not imply 
that it is correct, and in the presence of independent tests 
which do reject the model, the overall result should be 
to reject the model. This is expressed mathematically as 
follows. 

The results for the two tests, for the two subsamples, 
imply a confidence level 

1-P = l~ PMM] Po.[l^mi,^] PsifM] P5[t^max] (6) 

for rejecting T-Ci- Note that although the fii^o) and i^rnax 
tests seem to be independent (e.g. Figs |[ |7^, this has not 
been strictly proven, so the combined confidence levels for 
rejecting Hi may be slightly overestimated. 

The final contours in confidence levels for rejecting 
Hi{^lo, Xq) are shown in Fig. ^. 

These show that the null hypothesis of the absence 
of the large scale structure peak is rejected at at least 
the 1 — P > 84% level for nearly all pairs (r2o,Ao) 
with 0<fJo;S0.6 and at the 1 - P > 98% level for 
0.1 ^ rio ~ 0-5- Moreover, for a band running from (JIq ~ 
0.15, Ao = -0.1) to (no « 0.4, Ao = 1), HiiQa^Xa) is re- 
jected at the 1 - F > 99.9% level. 

In other words, if the matter density of the Universe is 
low, then the possibility that there is no large scale struc- 
ture peak at ^lss = 130 Mpc in the quasar sample is 
rejected, and it is rejected to high significance for the most 
favoured values of the (fio, Ao) pair: a low density hyper- 
bolic model with (Oq « 0.2, Ao = 0), or a low density flat 
model with {flo « 0.3, Ao = 1 — ilo)- 

For a flat, critical density model, {flo = 1, Ao = 0), Ti-i 
is not rejected. Could it be argued that if (ilo = 1, Ao = 0) 
is correct, then the rejection of Hi for low values of SIq is 
simply an artefact due to making a wrong assumption? 



Fig. 9. Confidence intervals for rejecting H2, the hypoth- 
esis that the large scale structure peak occurs at 1 /^Lgg , 
for the right ascension subsample. Shading styles are as 
for the previous figures, except that the confidence levels 
are two-sided, i.e. the four successively darker shadings 
are for 1 - P > 0%, 1 - P > 68%, 1 - P > 95% and 
1 - P > 99.7% respectively. 

A simple, quantified counterargument to this is the fol- 
lowing. If a peak really is present in the data, even though 
its amplitude may be low, then even for incorrect values 
of (r^o, Ao) it is likely that non-random frequencies can be 
detected. Since the search for the frequency starts just a 
little above 1 / ^Lgg ' then if a low value of fio is correct, the 
case of (57o — li ^0 = 0) may lead to the detection of a low 
frequency harmonic. This is indeed the case. By substitut- 
ing the word 'lower' for 'higher' in (iii) of the definition of 



Hi (Sect. 2.3), and recalculating the equivalent of Fig. 
the confidence level for rejecting Hi{^lo = l,Ao — 0) is 
found to be at the 1 - P > 98% level. 



3.2. 7^2 ■' What pairs (rio,Ao) are consistent with the 
frequency of the peak being at 1/ ^Lgg ^ 

Given that Hi is strongly rejected for interesting pairs 
of (f2o:Ao) values, what are the pairs (f2oiAo) which are 
consistent with ^Lgg ~ 130 ± 10 h^^ Mpc? 

For consistency, this question is formally answered by 
trying to reject the null hypothesis H2{^o, Xo), according 
to which the frequency of the peak is assumed to be at 
^LSS — 130±10/i~^ Mpc. Two-sided confidence intervals 
are used [a factor of two is inserted in front of the integral 
in Eq. (||)] , since the question of interest is now how close 
the best estimate of the frequency is to the hypothesised 
frequency, and both low and high frequencies would reject 
the hypothesis. 

The region in the (fJoiAo) plane for which H2 cannot 
be significantly rejected provides an estimate of the values 
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Fig. 10. Confidence intervals for 7^2, for the declination 
subsample, shading as for Fig. ^ 




0.5 



Fig. 11. Confidence intervals for TL2, combining the infor- 
mation from the two subsamples. Shading is as for Fig. ^. 

of r^o and Aq. In the present for most other meth- 

ods which provide significant constraints on the curvature 
parameters, there is a degeneracy between the latter, so 
that an estimate can only be provided for the relation 
between them, rather than for both parameters indepen- 
dently. A linear fit is used to describe this relation. 

Figs |9| and |l^ show the confidence levels for rejecting 
7i2(f^o,Ao) for the right ascension and declination sub- 
samples respectively. Neither subsample is sufficient on its 
own to provide a precise estimate of an (f2o: -^0) relation, 
although both reject to >,70% significance. 

However, the fact that they are independent subsam- 
ples implies that they can be combined via 

l-P = l- Po. [J^max] Ps [t^max] (7) 



to give Fig. ^ which shows that most values of (fJo, Aq) 
can be rejected to > 68% confidence if jfio — 0.25| >, 0.15. 

More precisely, a linear relation between the two cur- 
vature parameters can be fitted to the points in Fig. pi] for 
which the probability P of obtaining the observations is 
highest, using the 68% confidence limits as uncertainties, 
by linear regression of ilo function of Aq. 

This relation is 

r^o = (0.24 ± 0.04) -I- (0.10 ± 0.08) Aq. (8) 

The uncertainties here relate to the fitting procedure. 
Since the 1 a (68% confidence) uncertainties on for each 
value of Aq are not mutually independent, the uncertain- 
ties in Eq. ^ do not represent measurement uncertainties. 

Since the measurement uncertainty (including Ai^gg) 
is already expressed in the contours in Fig. |ll|, this is 
restored to the zero-point of the relation, giving 

VLo = (0.24 ± 0.15) + (0.10 ± 0.08) Aq, (9) 

where the zero-point uncertainty includes known measure- 
ment uncertainties and the slope uncertainty relates to the 
fitting procedure. 

Note that, as revealed by this relation, the point least 
rejected by the data for Aq = is much closer to the large 
ilo 68% confidence limit than to the low JIq 68% confidence 
limit, so the representation by Gaussian uncertainties is 
not an optimal approximation. However, if one deduces 
ilo = 0.24 ±0.15 for the case Aq = and assumes that this 
is a Gaussian 1 a uncertainty, then this will be sufficient 
for most applications, where one prefers overestimates of 
uncertainties to underestimates. 

Alternatively, the Aq = result can be written as fio = 

For a flat universe, i.e. Aq = 1 — fio; use of the un- 
certainties in Eq. (^) as Gaussian uncertainties and com- 
bination in quadrature yields fio = 0.30 ± 0.15, which is 
consistent with Fig. |l^. 

A large part of the uncertainty here is due to the boot- 
straps. For example, if the intrinsic measurement uncer- 
tainty due to the bootstraps is removed, and only the un- 
certainty AL^gg = 10 Mpc is used for the equivalent 
of Fig. |ll|, then the regions not rejected by the confidence 
level contours shrink considerably, and Aq < 0.4 is rejected 
at the 1 - P > 95% level. 

However, although bootstrap estimates of uncertain- 
ties provide an upper estimate to uncertainties, i.e. the 
true uncertainties may be smaller, it is prudent to retain 
the bootstrap estimate. 

3.3. Amplitude of large scale structure peak 

Since the main aim of the present study is to use the ^Lgg 
scale as a ruler for measuring the curvature parameters, a 
statistically robust estimate of the amplitude of the power 
spectrum peak used as a standard ruler is beyond the 
scope of this paper. 
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Indeed, as can be seen in Fig. ^ (upper plot), the am- 
plitude of the peak in the right ascension sub-sample is 
insufhcient (on its own) to significantly rule out the hy- 
pothesis of no peak at all, with respect to the Poisson 
simulations. Of course, the search for a local maximum 
does find that this is closer to vq than expected randomly, 
for a certain band in the {^o, Aq) plane. 

However, as a guide to what might be expected in 
future less sparse surveys, the amplitude in the declina- 
tion subsample, which has a stronger signal than that of 
the right ascension subsample, may be useful to quantify, 
though caution is recommended in the interpretation of 
this estimate. 

For the pair {Hq — 0.3, Aq — 0.7) (see Sect. a crude 
estimate of the amplitude expressed as a signal-to-noise 
ratio, 



A 



/(i^olobsvn) - (/(t^o|Poisson))_^ cr[/(t/o|Poisson)] ^_^^^ 



(/(i^o I Poisson)) 



(/(i^olPoisson)) 



is A« 1.7 ±0.5. 

This value is l ower than the corresponding value in 
Broadhurst et al. 's (1990) one-dimensional survey, for 
which A « 7 at the -Z^^gg scale (from fig. 2b of that pa- 
per), but is similar in order of magnitude to the density 
contrast values o f Sp/p>2.5 found in the Las C ampanas 
Redshift Survey ( Tucker, Lin fc Shcctman 199§| ). 

This is not surprising. Even for {Qq ~ 0.3, Aq — 0.7), 
a somewhat lower amplitude of density contrast can be 
expected at z ~ 2 relative to z ^ 0, though this could be 
compensated for (or under- or over-compensated for) by 
positive biasing of the quasar distribution relative to that 
of galaxies, if quasars turn on at the densest points where 
galaxies are most likely to interact and/or merge. 

A full analysis of the amplitude of the signal in analy- 
ses following the present one should potentially provide a 
useful constraint on models of quasar onset and lifetimes. 

3.4- Selection Effects 

The results above are strikingly consistent with the most 
recent expectations from independent observations re- 
garding large scale structure and the curvature param- 
eters: the power spectrum peak is present at ^lSS ~ 
130 ± 10 Mpc for the popular curvature pair {^Iq — 
0.3 ± 0.15, Ao==l-r2o). 

Could this just be a coincidence due to selection ef- 
fects? In the redshift direction, selection effects have long 
led to surprising results, though not to expected results. 



from quasar catalogues (see Scott 1991 and references 
therein). 

The angular scale corresponding to z ^ 2 and (fio ^ 
0.3, Ao = 1— ^lo) is ~ 2°. This is half an order of magnitude 
smaller than the size of a UK Schmidt plate (~ 6°). This 
is sufficiently small that the large scale structure scale is 
clearly smaller than the size of the plates, but not so much 
smaller that more subtle effects related to the plate size 



can be trivially excluded from contributing to the result 
found above. 

Possible angular selection effects, instrumental and/or 
astrophysical, in objective prism quasar surveys include 

(i) effects due to human subjectivity of selecting 'quasar- 
like' objects from the photographic plates 

(ii) not finding quasars in regions where the sky back- 
ground is noisy or there is no signal at all due to step 
wedges, large, bright galaxies or to bright stars and 
their ghosts (due to reflection from the secondary mir- 
ror support structure) 

(iii) not finding quasars close to other quasars/stars due 
to overlapping spectra 

(iv) differential apparent magnitude limits due to the vi- 
gnetting function of the telescope plus instrument ge- 
ometry 

(v) differences in apparent magnitude limits between 
plates 

(vi) differential apparent magnitude limits due to inter- 
vening dust 

(vii) mistaking quasars for stars in low projected number 
density open (star) clusters which happen to lie in the 
survey region, thereby missing quasars in those regions 

(viii) mistaking stars for quasars, which possibly explains 
the excess numbers of objects at some specific redshifts 
in Fig. ||. 

Problem (i) is avoided by application of a computer 
algorithm to digital scans of the Schmidt plates. The va- 
lidity of the precise quantification and relative weighting 



of the various 'quasar-like' criteria chosen by lovino et al 
(1996D to detect quasars could be debated, but since they 



are calculated automatically over the entire scanned re- 
gions, this aspect of human subjectivity is applied in a 
consistent, objective fashion across the plates. 

Problems (ii), (iii) to some extent (quasar-star over- 
laps), (vi), (vii) and (viii) are likely to be minimised by 
the choice of a very high galactic latitude region, i.e. the 
SGP. 

Problem (v) is corrected for in the control simulations 
by Poisson distributing points independently within the 
boundaries of ind ivid ual plates in each of the two sub- 
samples (see Sect. 



2.2 



The largest obvious contaminants listed in item (ii) 
which ocur in the fields studied here are a moderately 
bright star near a ~ l'',<5 ~ —30° (field 411, see ta- 
ble 1 of lovino et al. 1996| ) and a bright galaxy near 
a ~ 0''50™,(5 ~ 39° (field 295), which occupy less than 
about 0.1 sq. deg. and 0.2 sq. deg. respectively. Neither 
corresponds to what visually appears to be a void in Fig. |[ 

Counting the right ascension and declination anal- 
yses separately, the fraction of the solid angle bi- 
ased by these two objects is <, [(0.2) -I- (0.2 -|- 
0.1)] sq. deg./116.4 sq. deg. = 0.4% of the total solid an- 
gle included in the borders (as defined in Table [|) . This 
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Table 2. Angular separation of plate centres (table 1, 



lovino et al. 1996D from the four borders [a(l),(2), i5(l),(2) 
in table, from the borders in strictly increasing numerical 
a OT d values] and maximum centre-corner angle (3) of 
the fields for the right ascension (a) and declination {5) 
subsamples, defined here in Table |l|, in great circle degrees. 
Field numbers (#) are ESO/SERC field numbers. 



# : 295 296 297 295 351 411 

a(l) 2.37 2.55 2.53 2.37 1.73 1.40 

a(2) 2.43 2.45 2.47 1.09 1.97 2.51 

S(l) 2.27 2.26 2.25 2.27 2.77 2.77 



5(2) 2.23 2.24 2.25 2.23 2.23 1.73 



(3) 3.33 3.41 3.38 3.29 3.40 3.74 



is unlikely to be sufficient to mimic a large scale structure 
signal at ~ 2° in the full data set. 

Problem (iii) for quasar-quasar overlap could in prin- 
ciple cause a weakening of the real signal and not a false 
signal, although the low number density of the quasars 
implies this effect should probably be small. This is be- 
cause overlaps are expected to occur at around 70" in the 
dispersion direction and 3" in the orthogonal direction. If 
close quasar pairs are missed due to this effect, then the 
angular distribution measured will be smoother than the 
intrinsic angular distribution, i.e. filamentary type struc- 
ture would be less easy to detect than it should be. 

Quasar-star overlap, and problems (vi), (vii) and (viii) 
should be uncorrelated with the intrinsic quasar distri- 
bution and are more likely to weaken any genuine signal 
rather than mimic an expected signal. 

3.4.1. Differential magnitude limits? 

Problem (iv) could, in principle, provide the largest sys- 
tematic error. Fig. 3 of Dawe (1984 ) shows that from « 2° 
to ~ 4° from the centre of the UK Schmidt Telescope field, 
the apparent magnitude limit can become less faint from 
~ 0.02 mag to ~ 0.2 mag respectively. That is, less quasars 
would be detectable by a fixed search algorithm towards 
the edges and corners of the plate than in the middle. 

Table || lists the angular distances, in great circle de- 
grees, from the boundaries and the corners of the plate 
images, assuming that the plate centres listed in table 1 
of [lovino et al. (1996D are the correct (B1950) centres of 
the actual fields observed by those authors. 

The borders of the fields, as analysed here (Table |l|), 
are mostly about 2 — 2.5° from the centres, so the magni- 
tude limits should vary much less than 0.1 mag over most 
of the plates. The furthest corners from the centres are 
mostly at about 3.3 — 3.4° from the centres, so the mag- 
nitude limit should be about 0.06 mag brighter at these 
corners. 



Are differences in the (solid angular) quasar number 
densities due to these magnitude limit variations visible 
in Fig. ^? Possible voids in three of the —42° corners of 
fields in this figure, and in the (a = 0''42™,(5 = -37.5°) 
corner of field #351 (the middle declination field) could 
conceivably be related to the variation in magnitude limit. 

's quasars 



However, since the majority of lovino et al 



(see fig. 2 of lovino et al. 1996) are roughly uniformly 
distributed over an interval of about one magnitude in 
apparent magnitude, a change in the magnitude limit of 
0.06 mag could at most change the quasar number density 
by - 6%. 



Moreover, a small fraction of the quasars in the lovinc 



et al. sample have z > 2.4 and are not studied here. These 
are likely to have the faintest apparent magnitudes, so 
~ 6% is an upper estimate to the expected reduction in 
number density at the furthest corners. 

The apparent voids visible to the eye are presumably 
seen as voids because the number density is at least an 
order of magnitude lower than average. This is not expli- 
cable by a magnitude limit varying by 0.06 mag. 

In addition, in at least some of the fields, the num- 
ber density of quasars appears to be lower in 'voids' near 
the centres of the plates, where the exposure ought to be 
deepest. 

Figure ^ suggests that the S ^ —30° field might 
appear to have a circularly symmetric central concen- 
tration of quasars, apart from a cluster/filament near 

(a - 0^42"", S 31°). This is partly because the right 

ascension centres of the declination fields are offset from 
one another, which is why narrow limits in right ascen- 
sion were chosen. This could affect an analysis using the 
boundaries adopted here if that analysis were carried out 
in right ascension in the 'declination' subsample. But that 
is not the case here: only the distribution in declination of 
the declination sample was studied, so the right ascension 
offsets cannot affect the results. 

Another reason why problem (iv) is unlikely to be sig- 
nificant in the present study is that, as Dawe (1984) points 
out, departures from circular symmetry of the magnitude 
limit are below the measurement limits (0.01 mag) of his 
empirical estimates, and are not expected theoretically. 

The variations in quasar density around the field cen- 
tres appear about as strong as the variations leading to 
voids towards field edges or corners. For example, the void 
centred around (a = 0''50'", 6 = -37°) in field #351 (the 
middle declination field) is adjacent to a quite dense region 
around (a = 0''55™,(5 = —37°), which is about equidis- 
tant from the field centre. This strong variation cannot be 
caused by the circularly symmetric vignetting of the UK 
Schmidt Telescope. 

Finally, shifting several of the borders defined in Ta- 
ble |l] by 0.5° typically modifies the slope and zeropoint of 
the (r^o, Aq) relation of Eq. (^ by only about ±0.01. Since 
magnitude limits should, statistically, have some effect at 
the corners, even if small, this shows that uncertainty due 
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to differential magnitude limits appears to be negligible 
relative to the basic result. 



3.4.2. Selection effects: summary 

The possible selection effects discussed above do not seem 
sufficient to provide a non-cosmological explanation for 
the present results. It remains that t here could be unusua l 
systematic effects not reported by [ovino et al. (199^ , 
e.g. dust on patches of the photographic emulsion, bright 
supernovae, very bright comet or satellite trails. However, 
the overall similarity in the results for the right ascension 
and declination samples makes it unlikely that these could 
explain the principal results. 



3.5. Could the signal he confined to just a small part of 
the data set? 

Although Fig. 1^ provides qualitative evidence that the sig- 
nal comes from structures spread through the full data 
set, it could be possible that this subjective judgment is 
wrong. Could it be the case that a single structure pro- 
vides most of the signal, and that very little signal comes 
from the rest of the full data set? 

One possible approach to investigating this is by sub- 
dividing the full data set into smaller 'independent' sub- 
sets. Given the sparsity of the data set, subdivision into 
smaller subsets is likely to increase the noise. The i^LSS 
scale has a primarily statistical meaning, and random vari- 
ation in the estimate of the value of ^lss (something like 
the standard error in the mean) is likely to be greater in 
the smaller samples than in the larger samples. 

However, in order to provide an illustrative answer to 
the question, a subdivision of the two subsamples (right 
ascension and declination) was carried out, into the four 
redshift intervals 1.8 < z < 1.95, 1.95 < z < 2.10, 
2.10 < z < 2.25 and 2.25 < z < 2.40. These provide 
radial intervals which are close to the ^lSS scale, at least 
for (fio = 0.3, Ao = 0.7). Analyses across the full angular 
scale were carried out for Tii and 7^2 as before, but divid- 
ing the data into eight independent subsamples instead of 
just two. 

The result is that the hypothesis 7ii(0o, Aq) is rejected 
more strongly than before, so that the equivalent of Fig. || 
is rejection 1 — P > 98% everywhere in the domain of the 
(ilo, Ao) plane studied in this paper, and 1 — P > 99.9% 
over most of this. 

This shows that the presence of signals of some sort 
is stronger in the smaller subsamples than in the larger 
ones. 

However, the length scales at which the signals occur 
in individual subsamples are noisier in the smaller sub- 
samples, so that if no allowances are made for this, then 
the final result for 7^2 , the hypothesis that there is a signal 
at the expected scale ^lss every subsample, is that it 




Fig. 12. Confidence intervals for 7^2 (cf Fig. ^l]), for split- 
ting up both of the two subsamples into four redshift sub- 
samples each, and increasing the uncertainty in the power 
spectrum peak length scale to AL^gg = 20/i~^Mpc 
(Sect. 3.5). Shading is as for Fig. ^, i.e. the small patch 



suggesting a hyperbolic universe model is a 1 — F > 95% 
('> 2a') confidence interval region, and the rest is for 
1 - P > 99.7%. 




Fig. 13. Confidence intervals for 7^2 (cf Fig. |l^), for the 
four redshift subsamples of the declination subsample, re- 
taining Ai^Lgg — 10 Mpc (Sect. 3.5). Shading is as for 
Fig. H, i.e. the lightest shaded region is for the 1 — P > 68% 
confidence interval. 



is also rejected, to 1 — P > 99.7%, for the full range of the 
(r^o, Ao) plane covered. 

In order to find consistent solutions, either the uncer- 
tainty AL^gg needs to be increased, in order to allow for 
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the increased uncertainty per subsample, or some of the 
subsamples have to be (a posteriori) dropped. 

Since the number of objects per subsample decreases 
by roughly a factor of four, a doubling of AL^gg to 
Ai^Lgg = 20 Mpc should be sufficient to allow for 
consistency with the analysis for the full redshift range. 
This results in a small region which provides a 1—P > 95% 
solution in a small patch close to the (fio — 0-3, Aq = 
0.7) point, but suggesting a hyperbolic universe model 
(Fig. 111). 

Alternatively, retaining Ai^gg = 10 h^^ Mpc and 
just considering the four redshift intervals of the decli- 
nation subsample leads to Fig. |l3|. This shows a band in 
the (f2o, Ao) plane within the 1 — P > 68% contour, which 
is clearly consistent with (though narrower than) the cor- 
responding contour for the declination sample considered 
as a single sample (Fig. |l0|) . 

It is therefore clear that the signal is present in most 
redshift interval subsets of the full data set, but that be- 
cause of the increased noise, in order to combine those 
signals into a consistent solution, a more sophisticated 
technique than that presented here would be necessary. 



4. Discussion and conclusions 

It has been shown in this paper that the use of the power 
spectrum peak corresponding to large scale structure as a 
standard ruler in the tangential distribution of a homoge- 
neous quasar survey at z ~ 2 provides a new and indepen- 
dent method of constraining the curvature parameters. 

More precisely, what appears to be the optimal choice 
of homogeneous tangential quasar surveys publicly avail- 
able was chosen and analysed. The quasar data analysed 
were the right ascension and declination subsamples (Ta- 
ble |l|) of the [ovino et al. (199^ ) survey as provided in the 
Veron-Cetty fc Veron (1998| ) quasar catalogue. 

The null hypothesis 7ii(ilo, Ao) according to which no 
peak can be detected at -L^gg — 130 Mpc in the one- 
dimensional Fourier transform of the tangential proper 
motion distance distributions, and according to which the 
best estimate of the frequency of any such peak is ran- 
dom, is rejected to high significance for most 'interesting' 
pairs of (Oo,Ao) in the range (0.0 < flo < 1.1,-0.1 < 
Ao < 1.0). The highest rejection (at a confidence level 
1 — P > 99.9%) is for a band in the (fig, Aq) plane running 
from - (0.15, -0.1) to ~ (0.4, 1.0). 

Inversion of the frequency condition (iii) of Tii (fio , Ao) , 
by replacing 'higher' by 'lower', and recalculating confi- 
dence levels results in a rejection of 7ii(Jlo = IjAq = 0) 
at 1 - P > 98%. 

Since Tii is rejected, this implies that a power spec- 
trum peak is present in the data. 

Can the power spectrum peak be seen as a 'periodic- 
ity' in the data? Fig. allows the reader to judge this 
subjectively for the pair {Qq = 0.3, Ao = 0.7). 



^0=0.30, Xo=0.70, a 
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Fig. 14. Tangential length distribution for {Hq = 
0.3, Ao = 0.7), for the right ascension (upper panel) 
and declination (lower panel) subsamples (in bins of 
12.8 Mpc). The null hypotheses test relate to a power 
spectrum peak at iLgg — 130 ± 10 Mpc. Differences 
in number density between the plate boundaries are not 
corrected for, these are corrected for in the simulations. 



Consistently with the points at {Qq = 0.3, Aq = 0.7) in 
the upper panels of Figs || and 0, a periodicity of significant 
amplitude should be hard to detect in the right ascension 
sample, but disccrnable in the declination sample. This is 
the case in Fig. p^ . 

Moreover, although the redshift direction is potentially 
plagued by selection effects, the wedge diagrams in Fig. || 
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show clearly that large scale structure can be seen in the 
z — 9 plane, where 6 = a cos 6 or = S a.s above. Voids 
at a scale of around ^Lgg ~ 130 Mpc and difficult to 
explain by selection effects are visible in this figure. 

A best estimate for the values of (fJojAo) consis- 
tent with the occurrence of a large scale structure peak 
in the tangential quasar distribution is found by try- 
ing to reject the null hypothesis H2{^o,^o) according to 
which the best estimate of the frequency of a peak below 
l/d± = 0.01 /i Mpc~^ in the Fourier transforms occurs 
at 1/(Llss ± AivLgg). Bootstraps from the observational 
data set are used to robustly provide an upper estimate 
to the observational uncertainty. 

Based on a linear fit to the resulting estimates in the 
(fio, Ao) plane, the relation f^o = (0.24 ± 0.15) + (0.10 ± 
0.08) Ao is considered to best summarise the constraints 
on the two curvature parameters. For zero cosmological 
constant or flat models, the curvature parameter estimates 
arc (no = 0.24^°;?^, Aq = 0) and (f^o = 0.30 ± 0.15, Aq = 
1 — r^oi ) roopoctivoly. 




Th eac estimates of the density parameter arc in rc 
markable agreement with analogous estimates obtained 



Fig. 15. Confidence intervals from combining Fig, pi] with 
the relation O.Silo — O.6A0 = —0.2 ± 0.1 from Pcrlmuttei 
et al. (1999| ). Shading is as for Fig. ||. 



from the kinematics of galaxy clusters (e.g. Carlberg 
Yee & Ellingson 1997| ), collapsing galaxy groups (Mamon 
1993 ]j_ as well as f rom the baryonic fraction in clu sters 



0.55 < An < 0.95, 



(11) 



([White et al. 1993|; |Mohr, Mathiesen fc Evrard 19991) and 



groups ( Hcnriksen fc Mamon 1994 ). 

The estimate of the cosmological constant under the 
assumption of a flat model is also remarkably close 
to those which now tend towards a flat, cosmological 
constant dominated universe from faint galaxy number 
counts, the galaxy depletion curve behind clusters, the su- 
pernovae type la method ( Fukugitaet al. 1990|; Fort ct al 



1997yC ;hiba fc Yoshii 1997| ; [Perlmutter et al. 1999| ; [Riess ct 
al. 199§p and cosmic microwave background methods. Re- 
sults from the latter are being updated rapidly at the mo- 
ment, and are dependent on numerous assumptions which 
do not enter either the present method or the supernovae 
type la method, so are not discussed here. 

The present method is very independent of the super- 
novae type la method: in choice of astrophysical object, 
in the redshift range and in the difference between using 
standard candles versus standard rulers. What is the re- 
sult of combining the two? 

The combination of the present result with the relation 
0.8^0 + 0.6 Ao = -0.2 ±0.1 from [Perlmutter et al. (i999| ), 



where the error is modelled as Gaussian, results in Fig. 15. 

Although the two relations are not quite orthogonal, 
they have different enough slopes that the uncertainty in 
both is considerably reduced, such that a nearly flat model 
is implied without using any cosmic microwave background 
information. 

A linear fit to the 68% confidence contour in Fig. 15, 
for ilo as a function of (Ao — 0.7), results in: 

no = (0.30 ±0.02) -I- (0.57 ±0.11)(Ao- 0.7), 



where as before, the uncertainties from linear regression 
relate to the fitting procedure and underestimate the true 
uncertainties shown in the figure. 

In this case, the maximum uncertainty in Qq for a given 
value of Ao in the range above is cr(r2o) = 0.11. Restoring 
this to the relation as before yields: 

flo = (0.30 ± 0.11) + (0.57 ± 0.11)(Ao - 0.7), 

0.55 < Ao < 0.95. (12) 

This combined result strongly supports the possibil- 
ity that the observable universe satisifies a nearly flat, 
perturbed Friedmann-Lemaitre-Robertson- Walker model, 
where 'nearly' is quantified as ±0.1 in the two dimen- 
sionless curvature parameters. It does not require prior 
assumptions on the cosmological constant, nor does it use 
cosmic microwave background (CMB) data. 

An independent confirmation that no systematic errors 
are present in the data studied here (or else an indepen- 
dent estimate of the systematic errors) would obviously 
be very desirable. 

If further observations confirm that large scale struc- 
ture can be used as a standard ruler, then surveys such 
as the 2dF (2 degree field) and SDSS (Sloan Digitial Sky 
Survey) quasar surveys should provide confirmation of the 
present results within the next few years. 

The former risks systematic error from the fact that 
the angular scale corresponding to ^lsS ~ l"^*-* ^ 
10 h-^ Mpc at 2 - 2 is 2°: this is the field size of the 2dF 
survey. Careful correction for this, use of lower redshift 
quasars, and a much higher signal-to-noise ratio would 
help in obtaining convincing results from the former. How- 
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ever, since the SDSS is digital, it should presumably avoid 
this problem. 

The present method would be ideally suited to a deep 
right ascension survey (i.e. for a narrow band in declina- 
tion but a long, e.g. 60°, band in right ascension), which 
avoids any risk of dependence of any angular scale near a 
few degrees but includes a fair sample of quasars at z ^ 2, 
sampling a few hundred quasars per great circle degree. 
A blind, deep, slit spectroscopic, right ascension survey, 
'blind' in the sense of taking spectra of all objects in a very 
narrow declination interval, would provide possibly one of 
the fastest ways of obtaining constraints in the (r2o,Ao) 
plane which should suffer little from possible selection ef- 
fects. This offers a valuable cosmological project for the 
planned liquid mercury LZT (Large Zenith Telescope). 

An other survey to which the present method could 
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